Single-electron heat diode 
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We introduce a new functional nanoscale device, a single-electron heat diode, consisting of two 
quantum dots or metallic islands coupled to electronic reservoirs by tunnel contacts. Electron 
transport through the system is forbidden but the capacitive coupling between the two dots allows 
electronic fluctuations to transmit heat between the reservoirs. When the reservoir temperatures 
are biased in the forward direction, heat flow is enabled by a four-step sequential tunneling cycle, 
while in the reverse-biased configuration this process is suppressed due to Coulomb blockade effects. 
In an optimal setup the leakage heat current in the reverse direction is only a few percent of the 
forward current. 

PACS numbers: 73.23.Hk, 85.35.Gv, 07.20.Mc 

Understanding the thermal transport properties of 
nanoscale systems offers insight into fundamental physics 
as well as opens up possibilities for new applications. 
In a solid-state environment heat is generally carried by 
phonons and electrons, the latter offering much more flex- 
ibility in manipulation and measurement. Therefore in 
recent years there has been a great deal of interest to 
study heat control with mesoscopic electronics, including 
devices for cooling, 1 thermoelectric power generation^ 
and thermal rectification^ These advances suggest that 
in the near future thermal currents could be manipulated 
with a similar level of versatility as electrons today. This 
entails the construction of complicated heattronic circuits 
out of elementary components. One basic building block 
in electronics is a diode, and its heattronic counterpart, 
a heat diode, is a two-terminal device that allows en- 
ergy to flow only in one direction. More precisely, with a 
temperature bias in the forward direction, a heat current 
J + flows between the terminals, while reversing the tem- 
peratures makes the current drop to J_. For a proper 
diode operation these two values should be separated by 
at least an order of magnitude, that is, the rectification 
ratio J+/J_ should be at least of the order of 10. In this 
Rapid Communication we introduce a new heat diode 
design where this requirement is easily fulfilled. That 
should be compared to recent experimental demonstra- 
tions of asymmetric heat fiow=^~— where the ratio has been 
much below 2. We also emphasize that although there 
is a considerable body of theoretical literature on ther- 
mal rectification (see, for example, Ref. and references 
therein), it is mostly concerned with simple model sys- 
tems with no concrete realizations available, while the 
present proposal can be straightforwardly fabricated and 
operated with current experimental technology. 

Our proposed device is a single-electron heat diode, 
consisting of two semiconductor quantum dots or two 
metallic islands (collectively called dots), coupled to 
different reservoirs (L and R) through tunnel barriers 
(Fig. []}. The coupling between the dots is purely ca- 
pacitive so that there is no electron transport through 
the system. However, with a temperature difference be- 
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FIG. 1. (color online), (a) Schematic of the heat diode, shown 
here with a forward bias. Electrons hop in and out of the 
dots and exchange heat through Coulomb interaction U. (b) 
Circuit diagram for the diode with capacitances and external 
voltages indicated. Shaded areas correspond to the dots. 



tween the reservoirs, thermal current is generated by elec- 
trons tunneling in and out of the dots and exchanging 
energy through Coulomb interaction. This energy trans- 
port mechanism has been previously studied in Ref. in 
the context of a mesoscopic heat engine. Here we ex- 
tend the idea to produce a rectifying mechanism for heat 
currents. Let us fix the notation by stipulating that in 
the forward direction the right reservoir is hot and the 
left one is cold, with temperatures Tjj and Tc, respec- 
tively. The energy required for an electron to tunnel from 
reservoir a into the adjoining dot is E an , where n is the 
occupation of the other dot. With two gate voltages the 
level structure of the double dot can be tuned to such a 
regime where the occupation of each dot may only be 
or 1, and where the energies Ej, n are of the same order as 
Tc and the energies En n are much larger, at least as large 
as Th (Fig. [2]). Now one quantum of heat, equal to the 
Coulomb interaction energy U = Eli — Et,q = Em—EnQ, 
can be transported from right to left with the four-step 
cycle depicted on the left side of Fig. [H^b). The time- 
reversed cycle carries heat in the other direction, and the 
total current is the sum of these two contributions. When 
the temperatures are reversed, electrons from the right 
reservoir, now at the low temperature Tc, are unable to 
tunnel into the adjoining dot because the available ther- 
mal energy is much smaller than the required energies 
En n . Therefore the heat-carrying cycle is exponentially 
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FIG. 2. (color online), (a) Energy-level diagram of the 
double-dot system, (b) Sequential tunneling cycles transport- 
ing a heat quantum U between the reservoirs. The processes 
denoted by the thick arrows become suppressed under reverse 
bias. 



suppressed and heat flow is blocked. This asymmetric 
Coulomb blockade configuration is the origin of the diode 
effect. The general idea of producing rectification by cou- 
pling the left and right reservoirs to different transitions 
of the central system has been previously discussed in 
Ref. |. 

The rest of the paper is organized as follows. First we 
present the device model in detail and calculate the heat 
current in the sequential tunneling approximation. Then 
we show quantitatively how the diode effect arises in the 
proposed device, considering two realizations, metallic 
islands with a continuous spectrum and quantum dots 
with discrete states, and conclude that their behavior is 
essentially identical. We analyze the conditions for opti- 
mal diode operation and show how to obtain the required 
level structure with a double-dot setup. Finally we con- 
sider the experimental demonstration of the diode effect. 

The system is modeled with 



H = 



E 

x=L,R 



(H a + Hua + Hxa) + He 



(1) 



where H a — ^2 i£a Ejcjcj is the Hamiltonian for the 
reservoirs, H Da = J^ieDa £ DiC f Di c Di for the dots, and 
H Ta = J2i £a , 3 eDa tiCjjjCi + h.c. for tunneling between 
reservoirs and dots. Hq contains the Coulombic charging 
and interaction energies. We study a parameter regime 
where the charge states of the double dot can be trun- 
cated to the four lowest levels which we call (0, 0), (1, 0), 
(0, 1), and (1, 1). The ground state is labeled (0,0) and 
(tiljUr) is the state with n a excitations in dot a. The 
particles occupying the dots in the excited states can be 
either electrons or holes, and the Coulomb energy U is 
positive or negative depending on whether the two exci- 
tations in the (1, 1) state have equal or opposite charges. 
Sanchez and Biittiker— use a similar energy level structure 
to produce a heat engine. However, their device requires 
energy filtering with single-level quantum dots, whereas 
our diode principle does not depend on the dot density of 
states and is therefore applicable also to metallic islands. 

In the sequential tunneling approximation heat current 
is calculated in terms of the rates Tan for tunneling be- 



tween reservoir a and the adjoining dot, with n being 
the occupation of the other dot and a = 1 for tunnel- 
ing into the dot, a = for tunneling into the reservoir. 
The probabilities P nL n R to be in state (nL,n R ) are ob- 
tained from a steady-state master equation, and the heat 
current through the system, equal to the net energy ex- 
tracted from the right reservoir, is then given by^ 
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where the two terms on the second line correspond to 
the two cycles in Fig. Hfb), the normalization is T 3 = 

Ea^En.^o.i^M^r^ +rS), and a top bar 
denotes the other possible value, for example, = 1, 
L = R. Fermi golden rule gives the tunneling rates as 
= 27rr Q F Q ((-l)*£ cm ), where T a = £ j£q \U\ 2 S{ei) 
is the tunneling strength for junction a, assumed to be 
energy independent. For metallic islands with a con- 
tinuous spectrum the effective reservoir occupation is 
F a (E) = Ev Da n a (E), where v Da = J^ieDa 5 i £ Di) is the 
island density of states and n a is the Bose function. For 
quantum dots with a single discrete level SDa we have 
F a (E) = f a (E), the Fermi function for reservoir a. Sub- 
stituting these in Eq. @ gives for a symmetric structure 
(T L =T R = T) 



J = A- 1 UT[n R (U)-n L (U)]. 



(3) 



The diode effect is contained in the asymmetric weight 
function A, given for quantum dots as 



A _ 2[l + n L (U)}[l + n R (U)} g 
/lo(1 — /z-i)/iio(l — /fli) 

and for metallic islands as 

[a + n a (U)][a + n a {U)] 



a=L,R 
n,o-=0,r 



E an [8 na + n an ][n - {-l) n f a s na 



(4) 



(5) 



fa(E, 



n a (E an ), and v D = v DL 



where f a 
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The performance of the diode can be assessed by ana- 
lyzing two key quantities, the heat current under forward 
bias, J+, and the ratio of forward and reverse currents, 
J_l_/J_. Forward bias is defined as the right reservoir 
having the high temperature Tjj and the left reservoir 
having the low temperature XcS under reverse bias the 
temperatures are reversed. We start the analysis by mak- 
ing simplifications to Eqs. ([3]) and First note that 
the constant term in Eq. ((4]) is practically always small 
compared to the total A and can be neglected. The four 
Fermi factors in the denominator of A are then divided 
into two groups: those with a value larger than | are 
taken to be constants while the rest are approximated by 
exponentials. The device has two qualitatively different 
operating regimes, one with a positive and the other with 
a negative U, and we first consider U > 0. In this case 
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we have 1 — f a i ~ 1, giving the dependence of forward 



current on the energies E a0 as J + oc e 
rectification ratio as 



J+/J- 



,( T c 1 ~ t h 1 )( e ro-E lo ) 



t h and the 



(6) 



We see that {Er — Elo)/Tc should be large for efficient 
rectification while Elo/Tc and Ero/Th should not be 
much above unity for a large forward current. Also, the 
difference in the temperatures should preferably be in the 
range Th/Tc > 2. Since by definition (0, 0) is the ground 
state, both energies E a0 must be positive, and thus the 
optimal value of Elo is zero. The energy Ero can then 
be chosen to achieve the desired balance between forward 
current and diode efficiency. The interaction energy U 
has only a small effect on rectification but the forward 
current depends on it strongly. When U 3> Th, the Bose 
functions in Eq. © suppress the current exponentially, 
while in the limit U <C Th the Bose functions in Eq. §4§ 
give a similar suppression. Thus the optimal value of U 
is of the order of Th- 

When considering the other operating regime, U < 0, 
we can restrict attention to the interval Elo < \U\ < Ero 
since a large \U\ would push the (1,1) state below the 
ground state, while the above approximations for posi- 
tive U are actually valid for U > —Elo an d therefore the 
(uninteresting) case of \U\ < Elo is contained in the pre- 
vious discussion. Now we can approximate 1 — fm ~ 1 
while the three other Fermi factors in A are taken as 
exponentials. Then the current depends on the denomi- 

nator of A as J+ oc e T c t h ; and the rectification ratio 
is 



J+/J- 



= Per 1 — Th 1 )(E r0 — \U\) 



(7) 



Here \U\ has replaced Elo in both expressions and there- 
fore the device operation is almost independent of Elo 
in this regime. Numerical calculations show that the for- 
ward current is maximized when Elo ~ 2Tc. 

Above we considered only the quantum dot setup, 
characterized by Eq. However, the same approxi- 

mations and conclusions can be derived from Eq. ([5]) by 
taking into account only the dominant two terms of the 
sum, that is, those terms where both occupation factors 
in the denominator are close to zero. The conclusions 
can also be verified by a numerical analysis of the full 
equations, as is done in Fig. [3] We see that both positive 
and negative U can produce similar levels of rectification 
but the forward current is larger for U < and there- 
fore we concentrate on that regime. To have an efficient 
diode with J + / J_ > 10, we see that unless a tempera- 
ture bias of Th/Tc > 3 is available, the forward current 
cannot be larger than about 10 _1 T*Tc, where r* = T 
for quantum dots and = TvdEro for metal islands. 
On the other hand, if a current level of 10~ 2 T*Tc is suf- 
ficient, considerable rectification takes place already for 
Th/Tc < 2. 

The validity of the above analysis requires that the 
dots are not coupled too strongly to the leads. Since 
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FIG. 3. (color online), (a) Forward current J+ and the rectifi- 
cation ratio J+/J_ as a function of U and Ero. Black curves 
correspond to quantum dots (r» = T), red for metal islands 
(r» = TvdEro). Dotted lines are the estimates from Eqs. ([6]) 
and 0. The parameters are T H = 2.5 Tc, E L0 = 2Tb, 
Ero = 7.5 Tc (left panels) and U = -4Tb (right panels), (b) 
Maximum rectification as a function of thermal bias Th/Tc 
for fixed values of forward current. The three curves corre- 
spond to, from top to bottom, J + /T*Tc = 0.01 (black), 0.05 
(blue), 0.1 (red). For each data point the parameters Elo, 
Ero and U are optimized to give maximum rectification for 
the given current. Optimal values of Ero are shown in the 
inset, for Elo/Tc they are in the range 2 . . . 3, for U /Tc in 
the range —5.5 • • • — 3.5. (c) Charging energies for the double 
dot as a function of the right gate charge Qgr- The curves 
correspond to different charge states, with the relevant ones 
marked as [Nl,Nr], where N a is the number of electrons in 
dot a. Dotted lines highlight two values for Qgr where the 
heat diode could be successfully operated, that is, the left and 
right transitions are well separated, and the lowest four levels 
are well separated from higher levels. Here the left gate charge 
is Qgl = .48 e, the capacitance constants are </l = gR — .1, 
and the unit of energy is e 2 /(2C(l — <7l<?h))- 



we require that only the lowest four levels are accessible, 
the tunneling rates should be smaller than the energies 
of the neglected higher-lying states. The energy scale 
for the forbidden transitions is set by Ero & n d thus we 
must have 27rr* -C Erq. This coincides with the re- 
quirement that the junction conductances must be well 
below the conductance quantum. On the other hand, 
coupling strength must also be limited to prevent leak- 
age currents due to cotunneling effects. As already dis- 
cussed, the suppression of current under reverse bias is 
due to the fact that the low-energy particles are unable 
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to tunnel into the right dot. This blockade can be lifted 
by coherent two-electron processes where particles tun- 
nel through both junctions simultaneously and most of 
the required energy is contributed by the left reservoir. 
There are two such cotunneling processes: if the system 
is initially in the (0, 0) state, particles tunnel into both 
dots and the final state is (1,1); the other possible process 
leads from (1, 0) to (0, 1). In the supplementary material 9 
we present a detailed calculation of these processes and 
arrive at the conclusion that they can be neglected when 
r* is bounded as above. 

Let us now consider how to actually realize the re- 
quired energy level structure in a double-dot setup. For 
a metallic island, the tunneling energy E an is simply the 
difference of electrostatic charging energies of the initial 
and final states, for the quantum dot system one must 
also add the bare energy of the discrete single-particle 
level, EDa- Using the notation of Fig. 1(b), the charging 
energy for dots with charges Ql and Qr is 

E c h{Qh, Qr) = 

9l{Ql + Qgl? + 9r{Qr + Qgr? + 2g L g R Q L Q R (8) 
2C(1 - g L g R ) 

with interaction constants g a — C/(C + C a + Gc a ) and 
gate charges Q Ga = C Ga (V Ga - V a ) + g & C & {V a - V a ) + 
9aPG&^yGa — V a ). We see that the level structure is con- 
trolled by two gate charges which in turn depend on three 
voltages, and therefore in a practical implementation one 
of the voltages can be discarded. If in the ground state 
(0, 0) the dots have charges Qlo and Qro, then in the ex- 
cited dots the charges are Qlo + Ql and Qrq + qr, where 
the q a are — e or +e depending on whether the added 
excitation is an electron or a hole. From Eq. ([8]) we then 
see that the magnitude of U is determined by the capac- 
itances as \U\ = e 2 /[C(g2 9r — 1)] while the energies 
El$ and Er as well as the sign of U can be chosen with 
the gate charges. The energies E a0 are 

Elo = E c h(QLo + Ql, Qro) — E c h(QLo, Qro) — —£dl 

e 

Ero = E c h(QLo, Qro + ?r) — E^fQ^o, Qro) — —£dr 

e 



where the single-particle levels EDa, vanish for metallic is- 
lands. If the truncation of the system to the four lowest 
charge states is to be valid, all higher states must have 
an energy much above all E an . Of course the levels E an 
must themselves have the structure discussed above. Fig- 
ure EJc) shows two concrete examples of how to obtain 
appropriate level diagrams based on the above consider- 
ations. 

Experimental demonstration of the diode effect must 
take place at around 1 K or below. First, the requirement 
of large absolute bias, Tjj/Tc ~ 2, can only be attained 
at low temperatures. Second, since the device is based 
on Coulomb blockade physics, the charging energies at- 
tainable in mesoscopic structures set the maximum op- 
erating temperature to about 1 K. Third, in order to ob- 
serve the thermal current through the diode, heat trans- 
port by phonons must be suppressed, and this generally 
takes place at sub-Kelvin temperatures.— In a realistic 
setup we could have Tq = 100 mK, Th = 250 mK, and 
r* = 0.5 T c , giving r*T c « 10 fW. Referring to Fig. S 
we see that it is then possible to have a forward current 
of the order of 1 fW with the reverse current being a few 
percent of this value. The forward current is large enough 
to be measured with state-of-the-art thermometry*^ 

In summary, we propose a new device, a single-electron 
heat diode, which can be realized by a double-dot system 
connecting two electronic reservoirs at different temper- 
atures. The device and required operation scheme are 
routinely realized with currently existing technology. We 
explored the rectification performance of the device in 
detail and showed that with experimentally measurable 
current levels it is possible to have a rectification ratio 
well above 10, even up to about 100, making the present 
device the first concrete proposal for an efficient heat 
diode. 
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One of the authors (T.O.) acknowledges the Academy of 
Finland for financial support. 
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SUPPLEMENTARY MATERIAL: VIRTUAL TWO-PARTICLE TUNNELING RATES FOR THE 

SINGLE-ELECTRON HEAT DIODE 



As explained in the main text, heat flow under reverse thermal bias is blocked because the golden-rule rates for 
hopping into the right dot are exponentially suppressed. This blockade can be lifted by coherent tunneling of electrons 
through both junctions simultaneously, with most of the energy contributed by the hot left reservoir. There are two 
different processes which can lift the blockade: if the system is initially in the (0, 0) state, particles can tunnel into 
both dots and the system ends up in the (1, 1) state. If, on the other hand, we start with the (1, 0) state, one particle 
hops out of the left dot and another one into the right dot, leaving the system in the (0, 1) state. Here we calculate 
the rates for these virtual transitions and show that in the parameter range relevant for diode operation they are 
negligible compared to sequential processes. The rates for two-particle transitions can be calculated from the general 



expression 



r (2) 



2nJ2\(h\H T 

fb,ia. 



■H T \i a )\ 2 W la S(E fb -E la ), 



(1) 



where Ht = Htl + Htr is the tunneling Hamiltonian, Hq = H — Ht the uncoupled Hamiltonian, and Wi a the 
probability weight for the initial states. For the process (0,0) — > (1,1) the final state is related to the initial state 
by |/b) = c DL CLC DR c R\^ a ) ■ Note that this only applies to the case when the excitations for both dots are electrons. 
If, for instance, the excitation tunneling into the right dot is a hole, one should make the change c^ DR c R — > cdrc r . 
However, the formulas below will be identical in all the different cases. Writing the rate explicitly for metallic islands 
yields 

(2) I 

Tool,.!! = 2-kY l T r v dl v dr I de L de DL dE R de DR f L (e L ){l - /l(£dl)]/,r(£,r)[1 - /fl(£Dfl)]x 

(2) 

S(el + £r — £dl — £dr — Elo — E R0 — U). 



£l — £dl — Elo — ^Tio £r — £dr. — E 
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The two energy denominators correspond to the two different orderings in which the electron hoppings can occur. To 
regularize the divergences of the denominators we have introduced the inverse lifetimes F nLUR of the state (til, n R ), 
calculated with the lowest-order golden rule. Thus we have, for example, Tiq = r^jj + r^j*. In the case of quantum 



dots the corresponding rate is 
T^ 11 =2ttT l T r J de L de R f L (e L )f R (e R ) 



+ 



£l — Elq — f Tio £r — E R o — ^Toi 
For the other relevant cotunneling process, from (1,0) to (0, 1), the metallic system has 



S(e L +e R - E L0 - E m - U). 



(3) 



.(2) 



2ttT l T r u dl v dr J de L de DL de R de DR [1 - /i(Ez,)]/z,(eDi)/fl(efl)[l - Jr{£dr)}' 



1 



1 



£dl — £l — El 



£r 



--DR 



e R i — -irn 



(4) 



S{— £l + £r + £dl — £dr. + Elo — E R q), 



and analogously for the quantum dot setup. 

Equations @ and (HJ) suggest that for the metallic system it is useful to define new integration variables e' a — 
s a - £Da- Then with the help of the identities / Q (e' Q + e Da )[l - fai^Da)] = n a (e D a)[fa(£Da) - fa{^' a + £Da)] and 
f deoa [fa(^Da) — fa(s' a + ^Da)] = s' a the variables eo a can be integrated out. For both metallic and quantum dot 
cases the delta function can be used to eliminate another of the remaining integration variables, and then all the 
tunneling rates can be expressed with a single equation: 



r (2) _ 
1 A ~ 



2nT L T R U 2 J de F L {A ~ e)F R (e) 



1 



£ — E R o — ^Tq 



1 



E R i — ^Ti 



(5) 



where F a (e) — f a (e) for quantum dots and F a (e) = evo a n a {£) for a metallic system. For the (0, 0) — > (1, 1) transition 
the parameters are A = Elo + E R0 + U , T = r i, and r x = r 10 , and for the (1,0) — > (0, 1) transition they are 
A = E R0 - Elo, Tq = r 00 , and Tt = r n . 
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It is important to notice that Eq. ([5]) gives the total transition rates, including both sequential tunneling and 
cotunncling. 2 Sequential contribution comes from the resonance peaks at e = Er and e — E RQ , while the rest, due 
to the peak in the function Fl(A — s)Fr(s), is the cotunneling contribution. The sequential rate can therefore be 
calculated by taking the limit \e-E- §r|~ 2 ->■ 2*<J(e - E). We get 

r (i) r (i) r (i) r (i) r (i) r (o) r (o) r (i) 

r (2) _ 1 RQ L LI , LP Rl . r (2) _ 1 R1 L LI , LP RP fP \ 

1 00^11. sen — r > -p ! 1 10->01, sen -p T W 

i 01 1 10 1 1 



00 



This is exactly the result one would expect; for example, the first term in the (0,0) — > (1, 1) process is the rate of 
the (0,0) — > (0, 1) transition, multiplied by the probability that the next transition leads to (1, 1). The second term 
similarly gives the rate for the process (0, 0) — > (1, 0) — > (1, 1). 

To calculate the cotunneling rate we note that for exp(A/Tjf) > 1 we can approximate i*z,(A — e)Fn(e) w 
£ -a/t h ( e e(T c -t h ) _|_ e -e/T H -pi f Qr ^ ne q uan t um dot system. This function has a peak at |e| < Tc and it decays 
exponentially with 1 / (T^ 1 — Tjj ) and Th for positive and negative e, respectively. In this region we can approximate 
Eru — £ ~ -E^Kri an d therefore the cotunneling contribution can be extracted by setting \e — E — fT|~ 2 — > E~ 2 in 
Eq. ©. The result is 

r (2) _ ^T L T R U 2 T C a/Th 

for quantum dots. Similar considerations apply also in the metallic case, and the cotunneling rate is 

r (2) _ 2^Y L Y R v DL u DR U 2 T c A _ A/Tff 

i A,cot F 2 F 2 ~; n 2~Tc e ' 

(2) (2) 

The relevance of cotunneling can now be estimated by comparing the magnitudes of r^. ot and Tsoq- For a symmetric 
system, = T R = T, with U > Th, we have 

,(2) 



scq 



where T c g = T for quantum dots and T D s = -kYvdTc for metal islands. We have used the rates for the (1, 0) — > (0, 1) 
transition since the ratio would be smaller for the (0, 0) — > (1, 1) transition. For U < —Th the latter process dominates 
and we have 

Es°t *r cS U 2 Tc Elo/Th gHifT -i_ r -i^ 

rL 2 q ^o^i " ( j 

For the parameter values corresponding to the plots in Fig. 3 of the main text, with r c ff < Tc, these ratios are smaller 
than 10%, and therefore we conclude that cotunneling effects are negligible for the proposed device. 
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